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Abstract 



We calculate the collective charge density excitation dispersion and spec- 
tral weight in bilayer semiconductor structures including effects of interlayer 
tunneling. The out-of-phase plasmon mode (the "acoustic" plasmon) develops 
a long wavelength gap in the presence of tunneling with the gap being pro- 
portional to the square root (linear power) of the tunneling amplitude in the 
weak (strong) tunneling limit. The in-phase plasmon mode is qualitatively 
unaffected by tunneling. The predicted plasmon gap should be a useful tool 
for studying many-body effects. 
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Collective charge density excitations (or, equivalently plasmon modes) in bilayer struc- 
tures have attracted a great deal of theoretical and experimental attention over the last 
sixteen years ever since the existence of an undamped acoustic (i.e., with a long wavelength 
dispersion linear in wave vector) plasmon mode was predicted |]J in semiconductor double 
quantum well systems. In an uncoupled bilayer system, ignoring any interlayer Coulomb 
interaction, each layer can support a two dimensional (2D) plasma mode |2| with a long 
wavelength (q — > 0) dispersion uo ~ g 1//2 where q = |q|, and q is the 2D wave vector. When 
the two layers are near each other (separated by a distance d in the z-direction with the 2D 
layers in the x-y plane), the 2D plasmons are coupled by the interlayer Coulomb interaction 
leading to the formation |j] of in-phase and out-of-phase interlayer density fluctuation modes 
u)±: an out-of-phase acoustic plasmon mode, uo_ ~ q, where the densities in the two layers 
fluctuate out of phase with a linear wave vector dispersion, and an in-phase optical plasmon 
mode, uj + ~ g 1 / 2 , where the densities in the two layers fluctuate in phase with the usual 2D 
plasma dispersion. These u± modes have been observed p,|J in double layer semiconduc- 
tor quantum well systems via inelastic light scattering spectroscopic experiments, and the 
observation of the uo_ ~ q mode in GaAs-AlGaAs multilayer systems is, in fact, the only 
unambiguous direct experimental observation of an acoustic plasmon mode in solid state 
plasmas in spite of the theoretical literature on the subject going back more than forty years 

In this Letter we consider the experimentally relevant issue of the collective mode dis- 
persion in bilayer structures in the presence of significant interlayer quantum tunneling, all 
earlier work on the subject having dealt with the zero tunneling limit. It is well-known that 
tunneling introduces qualitatively new physics by introducing a new energy scale, the 
interlayer tunneling energy, in addition to the Coulomb energy and the intralayer kinetic 
energy. We find that tunneling significantly affects the out-of-phase mode, qualitatively 
modifying its long wavelength dispersion to co>_ ~ (A 2 + C\q + C^g 2 ) 1 / 2 where A defines 
the plasmon gap A = uo_(q = 0) which depends nontrivially on the 2D electron density 
n and the interlayer tunneling amplitude t. In particular we obtain the interesting result 
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that A ~ t or i 1 / 2 depending on whether the interlayer quantum tunneling is strong or 
weak respectively. We also find that, in contrast to the no-tunneling situation when the u + 
in-phase mode exhausts the plasmon spectral weight in the long wavelength limit (and the 
co>_ acoustic mode carries significant spectral weight only at finite wave vectors), the out-of- 
phase mode uo_ may carry significant spectral weight in the presence of tunneling even in 
the long wavelength limit and may be easily observable via inelastic light scattering spec- 
troscopy or frequency-domain far infrared (or microwave) spectroscopy J2]. We note 
the somewhat non-intuitive result that finite tunneling in fact converts the out-of-phase 
"acoustic" plasmon mode (in the t = situation) to an "optical" plasmon mode (in the 
t 7^ situation) by producing a finite plasmon gap u_(q = 0) = A whereas the original 
"optical" plasmon (for t = 0) in-phase mode u+ becomes the "acoustic" plasmon mode, 
albeit with a g 1//2 long wavelength dispersion, in the sense that U-(q = 0) vanishes in the 
presence of finite tunneling. The situation in the presence of tunneling is therefore similar 
to the familiar phonon terminology where the optical phonon (which has a finite energy at 
zero wave vector) corresponds to the out of phase intracell ionic dynamics and the acoustic 
phonon (with vanishing long wavelength energy) corresponds to the in-phase intracell ionic 
dynamics. 

The collective mode spectrum is given by the zero of the dynamical dielectric function of 
the system, which for a bilayer system in the presence of a finite interlayer tunneling ampli- 
tude t becomes a tensor eyj m of the fourth rank where i, j, I, m=l or 2 is the layer index with 
1,2 denoting the two layers. The dielectric function eijt m (q,u) = Su5j m — Viji m (q)Hi m (q, uS) 
is obtained within the mean field random phase approximation (RPA) in our theory where 
the deltas are Kronecker delta functions and Vij\ m is the intralayer/interlayer Coulomb in- 
teraction matrix element with Hi m as the irreducible noninteracting electron polarizability 
function. It is convenient to use one electron energy eigenstates E± (we take h — 1 through- 
out this paper) as the basis set rather than the layer index since the latter is not a good 
quantum number in the presence of tunneling. The energy levels E± = e(k) ± t, where 
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e(k) = k 2 /2m is the parabolic one electron 2D kinetic energy in each layer and t is the tunnel- 
ing strength, are the usual symmetric and antisymmetric states in the presence of tunneling 
with a single particle symmetric-antisymmetric (SAS) gap given by As^s = E + — E- = It. 
In the SAS representation the collective mode spectra become decoupled by virtue of the 
symmetric nature of our double quantum well system (i.e., both layers identical with equal 
electron density), and the collective density fluctuation spectra are given by the following 
two equations for the in-phase and the out-phase plasmon modes u± respectively: 

€+(<?, u) = 1 - v + (q) \n++(q, u) + II__(g, to)} = 0, (1) 

and 

e4q,uj) = 1 - v^q) [II + _(g,u;) + U-+(q,u)] = 0. (2) 

In Eqs. (|l]) and (0) the Coulomb interaction matrix elements are given by v±(q) = v\(q) ± 
V2(q) where Vi^{q) are respectively the intralayer and interlayer Coulomb interaction matrix 
elements. We use the simplest model for the Coulomb interaction (at no loss of generality) 
assuming the intralayer Coulomb interaction to be purely a 2D Coulomb interaction (and 
thus neglecting subband effects in each layer, which is entirely justified in most experimental 
situations where the intralayer intersubband energy is much larger than As^s) - subband 
effects can be trivially incorporated by using the appropriate subband form factors 0. 
For this simple model, V\(q) = 2ne 2 /{Kq)\ v^q) = i>i(?) exp(— qd), with k as the (high 
frequency) background lattice dielectric constant. Finally, H a p(q, u>) in Eqs. (|l|) and (|J), 
with (a, (3) = (+,—), are the noninteracting SAS polarizability functions within our RPA 
theory: 

tt ( n f d 2 k Uk + qj-fpjk) 

-2 J J^ w + Ea{k + q) _ E ^ k y ( 3 ) 

where f a p are Fermi occupancy factors (we restrict ourselves to T = OK in this paper), and 
the factor of two in the front arises from spin. 

Solving Eqs. (0)-(0) we obtain the collective density fluctuation spectra of the coupled 
bilayer system. In the absence of tunneling, t = 0, one has E + = EL = e(k), and one then 



recovers in a straightforward fashion the well-known optical (u + ~ q 1 / 2 ) and acoustic ~ 
g) plasmons of a bilayer system without any electron tunneling. It is, in fact, straightforward 
(but quite tedious) to obtain analytically [from Eqs. (TJ)-®] the long wavelength (q — > 0) 
plasma modes of the coupled bilayer system including effects of interlayer tunneling. We 
obtain in the long wavelength limit the following results: 

9 / \ 2vre 2 iV . . 

ul(q ^ = q, 4 

w_(g->0) = A 2 + dg + C 2 g 2 , (5) 

where N = 2n is the total 2D electron density (n being the electron density per layer), 
C*2 > 0, and 

A 2 = A SAS H ( n + - ^-)(gTF^)A S AS 

m 

C 1 = -(A 2 -A 2 SAS )± (6) 

where = 2me 2 /n is the 2D Thomas-Fermi wave vector and n± = n ± n c is the electron 
density in the symmetric/antisymmetric ± level. Here, n c = (m/2n)AsAS for n > n c when 
both symmetric and antisymmetric levels are occupied (i.e., the 2D Fermi energy Ep > 

as), and n c = n for n < n c when only the symmetric level is occupied (i.e., Ep < A$as)- 
Note that the in-phase mode uj + is unaffected in the long wavelength limit either by finite 
tunneling or by level occupancy, and depends only on the total 2D electron density N, 
following the standard 2D plasmon dispersion. The positive coefficient C*2(gTF, d, Ep, As^s) 
is not shown for brevity. 

The most important qualitative feature of the plasmon dispersion in the presence of 
interlayer tunneling is that the out-of-phase mode a>_, which is purely acoustic in the absence 
of tunneling {uj_ ~ q if A S as — 0), develops a plasmon gap at q = in the presence 
of nonzero tunneling. It is easy to see from Eqs. (§]) and (g) that this plasmon gap A, 
oj-(q = 0) = A, has the following behavior 



A ~ A SAS or J A SAS} (7) 



depending on whether the interlayer tunneling is strong (Asas/Ef 3> qrFd) or weak 
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(Asas/Ef <C qTFd). It should be emphasized that the strikingly non-intuitive A S ' AS depen- 
dence of the collective mode gap (on the square root of the single particle gap) is purely a 
Coulomb interaction effect, which dominates the collective excitation spectra in the weak 
tunneling situation. Finally, it is interesting to note that the first order dispersion correction 
to the out-of-phase plasmon gap is negative. 

In Figs. [I] - H] we present our numerical results for the collective excitation spectra of the 
coupled bilayer system without restricting to the long wavelength limit. In our calculations 
we have used RPA and also the so-called Hubbard approximation (HA) which includes a 
model static local field correction || to the noninteracting RPA irreducible polarizability 
[Eq. @]. In general, HA should be a better approximation than RPA at lower electron 
densities (and higher wave vectors) although, being an uncontrolled approximation, the 
quantitative improvement in HA over RPA is unknown. In Figs. [I] and |2| we show our 
calculated collective mode dispersions for different electron densities (with a fixed interlayer 
separation of d = 200A) in both RPA and HA, and for both zero and non-zero tunneling 
(Ag^s = 0, A S as = 1 meV). The important features to note in these results are: (1) In 
general both tunneling and local field correction have little effect on the in-phase mode u + , 
which even at low densities seems to be well described by the long wavelength RPA formula 
[Eq.(||)]. This is a direct consequence of the /-sum rule - the u + mode being the effective 2D 
plasma mode of the system is robust and insensitive to many body and/or tunneling effects. 
This somewhat disappointing result is however important because it states emphatically that 
all efforts to study many-body effects by studying the usual in-phase 2D plasma dispersion 
are doomed to failure, a fact already empirically known to experimentalists. (2) The out- 
of-phase uj_ mode is strongly influenced by tunneling and local field effects, and could in 
principle be a sensitive experimental tool in studying many-body effects |J. (3) Local field 
effects in general reduce the U- frequency, and at low densities this could even lead to a 
complete suppression of the u_ mode. In the absence of tunneling this complete suppression 
of uj- mode occurs below the critical density n c = q^ F /[8n(l + qrFd) 2 } in the HA. Tunneling, 
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however, has the dramatic effect of making the suppressed u- mode reappear when the layer 
separation d increases above d c = l/2mAsAs ~ ^-/qtf even if n < n c . (Note that, within the 
RPA, this Landau damping induced suppression of U- mode does not happen, and the U- 
mode exists for small wave vectors at all densities.) (4) Tunneling opens up a zero frequency 
plasmon gap in the u;_ mode. In RPA the calculated gap A = uo-(q = 0) is always greater 
than As as [Eq. ©]. In the HA, however, A could be above or below As as, an d is given by 



A 2 = AL5 + - n_){q TF d) (l - ^Tj) A SAS . (8) 



7T , „ A 1 

If the total density N < N c = 1/ (87rd 2 ) , the plasmon gap A is less than the single particle 
gap Asas, and we have the interesting situation of a collective plasmon mode being lower in 
energy than the corresponding single particle energy. (5) Tunneling leads to a weak negative 
dispersion of the u;_ mode at long wavelengths (i.e., G\ < 0). (6) Finally, a very interesting 
finding (Fig. |2|) is that at low densities the u;_ mode is, in fact, more stable in the HA (than 
in RPA) because it lies below the single particle continuum in the presence of tunneling. 
In the absence of tunneling the uj- mode may reappear at large wave vectors beyond the 
single particle continuum although it is completely suppressed (by Landau damping) at long 
wavelengths. This reappearance of the uo_ mode at large wave vectors (for As as = 0) is 
purely a local field effect, and has earlier been attributed Jl0| to a charge density wave 



instability in bilayer structures. Tunneling, according to our theory, stabilizes the long 
wavelength uj- mode which now lies between the two (symmetric and antisymmetric) single 
particle continua, and is in fact more stable in HA than in RPA. Tunneling, therefore, 
opposes the presumptive charge density wave instability |10[ . 



In Fig. |3] we show our calculated spectral weight or the dynamical structure factor 
S(q, q z , uj), which is given 0] by the imaginary part of the density- density correlation function 
where q z is the probe wave vector normal to the layers 0,^,0], for the collective modes in 
coupled bilayer systems both in RPA and HA for finite tunneling. For q z d = 0.0 only 
the in-phase mode uj + carries any weight, but for finite q z d the out-of-phase mode carries 
substantial spectral weight even at long wavelengths (q — > 0). The most important message 



here is that both u± modes in general carry finite spectral weights and should be observable 
in resonant inelastic light scattering [[3],|J and far infrared optical spectroscopies. 

We point out that the ui- mode, which presumably becomes the Goldstone mode in the 
symmetry-broken phase, could be used as an experimental probe to search for the theo- 
retically predicted |TTJ interlayer-spontaneous-phase-coherent quantum state in low density 
bilayer structures. It is interesting in this context to point out that it has been claimed [12| 
that the local field correction by itself within the so-called quasilocalized charge approxima- 
tion |12) , even in the absence of any interlayer tunneling, could lead to the opening of a long 



wavelength gap in the out of phase mode in bilayer systems. 

We note that the main approximations of our theory, namely treating the tunneling 
amplitude t as a phenomenological parameter and using the Hubbard approximation for the 
local field correction, should not in any way affect our qualitative conclusions. In particular, t 
decreases with increasing layer separation d in a calculable way, and therefore for large d one 
always recovers the no-tunneling limit of a vanishing plasma gap. While there is no general 
consensus on the best possible local field corrections in electron liquids, it is well-known || 
that the Hubbard approximation gives quantitatively similar results to more sophisticated 
local field corrections involving self-consistent approximations ||. We emphasize that the 
existence of a plasma gap in the presence of tunneling and its asymptotic dependence on t 
are independent of our approximations. 

In conclusion, we establish in this paper that bilayer semiconductor double quantum well 
structures should be a useful tool for studying the interplay between tunneling and many- 
body effects on collective mode dispersion - in particular, the dispersion of the out-of-phase 
collective mode in the presence of tunneling should show interesting observable many-body 
effects. 

This work is supported by the U.S.-ARO and the U.S.-ONR. One of us (S.D.S) thanks 
J. Eisenstein for asking an important question. 
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FIGURES 

FIG. 1. RPA plasmon dispersion for (a) n = lC^cm -2 , and (b) n = 10 10 cm -2 . Solid (dashed) 
lines represent plasmon dispersion in the absence (presence) of interlayer tunneling. Inset shows 
that the finite wave vector minimum of the out-of-phase mode uj- in the presence of tunneling. 
We use parameters corresponding to GaAs quantum wells: m = 0.067m e , k = 10.9, and d = 200 A, 
Asas = 1-0 meV for all our figures. Shaded region indicates the single particle excitation (SPE) 
Landau damping continuum. 

FIG. 2. Local field effects on the plasmon dispersion for density n = 1.0 x 10 9 cm~ 2 : (a) in 
the absence of tunneling, (b) in the presence of tunneling (As^s = 1.0 meV). Inset shows the 
reappearance of the out-of-phase mode w_ near q = 2fcp. 

FIG. 3. The dynamical structure factor for (a) n = 10 10 cm~ 2 and (b) n = 10 9 cm~ 2 in the RPA 
(solid lines) and HA (dashed lines) for finite tunneling (As as = 1-0 meV). Here, thin (thick) lines 
correspond to q z d = 0.0 (q z d = 7r/4). 
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